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The renormalized expectation value of the energy-momentum tensor for a scalar field with any 
mass m and curvature coupling £ is studied for an arbitrary homogeneous and isotropic physical 
initial state in de Sitter spacetime. We prove quite generally that {T a b) has a fixed point attractor 
behavior at late times, which depends only on m and £, for any fourth order adiabatic state that 
is infrared finite. Specifically, when m 2 + > 0, (T a b) approaches the Bunch-Davies de Sitter 
invariant value at late times, independently of the initial state. When m = £ = 0, it approaches 
instead the de Sitter invariant Allen-Folacci value. When m = and £ > we show that this state 
independent asymptotic value of the energy-momentum tensor is proportional to the conserved 
geometrical tensor ^H a b, which is related to the behavior of the quantum effective action of the 
scalar field under global Weyl rescaling. This relationship serves to generalize the definition of the 
trace anomaly in the infrared for massless, non-conformal fields. In the case m 2 + £R = 0, but m 
and £ separately different from zero, (T a b) grows linearly with cosmic time at late times. For most 
values of m 2 and £ in the tachyonic cases, m 2 + £R < 0, (T a b) grows exponentially at late cosmic 
times for all physically admissable initial states. 

I. INTRODUCTION 

Quantum field theory in curved spacetimes does not contain in itself a unique specification of the quantum state of 
the system [[[]. Even in Minkowski spacetime, where the existence of the Poincare group singles out a special state, the 
Minkowski vacuum, it is certainly of interest to consider states that are non-invariant under Poincare transformations, 
since they contain all the information about the physical excitations and dynamics of the theory. Such non-vacuum 
states are also necessary in a general initial value formulation of the back-reaction problem in both curved and flat 
spacetimes. In flat space the initial value problem for arbitrary physically allowable states has been formulated and 
studied for both QED and scalar $ 4 theory in the large N limit, principally for time varying but spatially homogeneous 
mean fields ||[|. 

The simplest situation in which the back-reaction problem can be studied in curved spacetime is that of a free scalar 
field in a spatially homogeneous and isotropic Robertson- Walker (RW) cosmology, where the geometry is characterized 
by just one non-trivial function of time. The wave equation for a free scalar field in such a geometry can be separated 
and expressed in terms of a complete set of time dependent mode functions. The general initial value problem is 
specified by giving initial data for this complete set at a given initial time. The back-reaction of the quantum scalar 
field(s) on the RW geometry can be studied then by constructing the renormalized expectation value of the energy- 
momentum tensor (T a {,) of the field(s) and solving (numerically) the semi-classical Einstein equations, augmented by 
higher derivative terms required by renormalization As in the flat space examples, this semi-classical back- 

reaction problem becomes exact in the large N limit, with N the number of identical scalar fields Jl0| . 

As a prelude to the dynamical back-reaction problem in cosmological spacetimes it is necessary to study non- vacuum 
states first in fixed RW backgrounds. The maximally symmetric de Sitter spacetime is of particular interest. Most 
previous work has focused on maximally symmetric 0(4, 1) de Sitter invariant states or the special 0(4) invariant 
state found by Allen |0|. Since the universe is not globally 0(4, 1) invariant, a more generic set of initial conditions, 
consistent only with RW symmetry and general principles of renormalization of {T ab ) is required for cosmology. The 
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investigation of these much weaker requirements and specification of the general initial value problem for back-reaction 
calculations was initiated in Ref. Jjjl , 

In this paper we study the behavior of the renormalized (Tab) for arbitrary physically admissable spatially homo- 
geneous and isotropic states in a fixed de Sitter background. We argue in Section III that such states must be fourth 
order adiabatic states |j| that also possess an infrared finite two-point function. In de Sitter space the wave equation 
for free scalar fields can be solved exactly for arbitrary values of the mass and the curvature coupling. Its solutions 
depend only on the wave number k of the mode and the parameter v 2 — j — m 2 a 2 — 12£, with R = 12a -2 the constant 
scalar curvature of de Sitter spacetime. For tft{v) < |, corresponding to m 2 + £R > 0, we prove that for all UV and 
IR physically allowed initial states the renormalized value of (T a b) at late times asymptotically approaches that of the 
Euclidean or Bunch-Davies de Sitter invariant state 13 - 16 1 . The conformally invariant scalar field (m = 0, £ = g) 
falls into this class. 

The case v = | corresponding to m 2 a 2 + 12£ = is more delicate. If m and £ are separately zero (the massless, 
minimally coupled case), then we prove that the renormalized (T a b) for all physically admissable states approaches the 
Allen- Folacci de Sitter invariant value [p7|-|l9|. Numerical evidence for this result was found previously in Ref. p2"l| . 
In this paper we provide an analytic proof that late time attractor behavior occurs for all physically admissable RW 
states, when m 2 > and £ > 0. If m 2 a 2 + 12£ = but m 2 and £ are not separately zero (so that one of them 
is negative) we prove that (T a b) grows linearly in RW comoving (cosmic) time without bound, and this asymptotic 
behavior is independent of the state of the field. Finally, and in contrast, in the case v > |, corresponding to 
m 2 +£R < 0, (Tab) depends sensitively on the state and, for most values of m and £, grows exponentially at late times 
for all states. This case is of considerably less physical relevance, since it corresponds to a tachyonic field theory with 
no stable vacuum state. 

The asymptotic approach of (T a b) to a de Sitter invariant form, independently of the lower symmetry of the 
initial data when ?ft(v) < | is a striking result. Certainly no such attractor behavior of (T a b), independent of initial 
conditions occurs in Minkowski space for any mass. One may regard this result as a kind of cosmic "no hair" theorem 
for scalar quantum fields in de Sitter space. For 9?(i/) < | it is in accord with one's classical intuition that any 
initial energy density satisfying the weak energy condition (e + p > 0) is redshifted away by the exponential de Sitter 
expansion, although as we will see, the redshifting of the quantum (T a b) is not that of classical matter or radiation. 
At asymptotically late times what is left behind is a kind of frozen quantum vacuum energy "condensate," satisfying 
the de Sitter invariant equation of state p = —e. This result justifies the choice of the Bunch-Davies vacuum in 
calculations of quantum fluctuations of free fields, i.e. without back-reaction, in a long-lived de Sitter expansion 
phase of inflationary cosmological models. For f=|,m = £ = 0, the approach of (T a b) to the de Sitter invariant 
Allen-Folacci value is perhaps more surprising. As shown in Section IV one expects the leading order contribution 
of the modes to (T a b) in this case to be constant in comoving time at late times. In fact this occurs if m and £ are 
not separately zero for all the modes. However when m and £ are both zero, the leading order contributions to (T a b) 
of all the modes except the spatially homogeneous one, for an arbitrary physically admissable state have exactly zero 
coefficient, the subleading contributions redshift away, and we are left only with the with the de Sitter invariant Allcn- 
Folacci constant value at late times. The finite difference from the Bunch-Davies value may be attributed entirely to 
the constant behavior of the spatially homogeneous mode contributing to the vacuum energy condensate in de Sitter 
space. 

In all those cases for which $t(v) < | when (T a b) approaches a de Sitter invariant value at late times, the quantum 
expectation value loses all its initial state dependence and hence its asymptotic value must be determined purely by 
the background geometry. When the mass of the field vanishes, the existence of only one covariantly conserved local 
geometrical tensor of adiabatic order four in de Sitter space, namely ^H a b given by Eq. (3.2) below, permits us 
to identify the asymptotic value of (T a b) in the vacuum energy condensate with this tensor. Since ^H a b cannot be 
derived by variation of a covariant local action, but corresponds instead to a certain well defined non-local term in 
the quantum effective action |pof , the asymptotic vacuum energy condensate of the quantum field is determined by 
the global or extreme infrared properties of de Sitter space. The form of the non-local effective action is determined 
by the trace anomaly in conformally flat spacetimes. Since the approach of (T a b) to a de Sitter invariant value occurs 
for all £ > 0, the existence of this term in the effective action for massless fields is much more general than the strict 
definition of the trace anomaly in the conformally invariant case. Hence the asymptotic late time behavior of (T a b) 
in de Sitter space can be used to define a generalized trace "anomaly" coefficient in the massless but non-conformally 
invariant cases, £ ^ \. As wc show by consideration of the covariant C function method flj, this coefficient is exactly 
the same as that which determines the infrared response of the vacuum condensate to global Weyl rescalings. Hence 
the significance of the state-independence of the vacuum energy condensate in de Sitter space is that it determines 
certain conformal properties of non-conformal field theories in the extreme infrared. 

The paper is organized as follows. In Section [H] we give the expectation value of the energy-momentum tensor as 
a mode sum for an arbitrary homogeneous and isotropic, physically admissable state with a non-zero initial particle 
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number of the scalar field. In Section III we analyze the late time behavior of expectation value of the energy- 
momentum tensor in de Sitter space in flat spatial sections and show that it appr oaches the de Sitter invariant 
Bunch-Davies value for all 5i(z/) < §, independently of the initial state. In Section IV we analyze the limit !ft(V) — > | 
in closed spatial sections in order to keep careful track of the spatially homogeneous mode in a discrete basis, and 
show how the Allen-Folacci fixed point at late times is obtained for the massless minimally coupled field. We also 
investigate the asymptotic behavior of the energy-momentum tensor for arbitrary mass and curvature coupling when 
v > |. In Section [y| we illustrate the analytic results with numerical studies, investigating in particular the interesting 
case when v is slightly smaller than | . We find that for many states when v is only slightly smaller than | the energy- 
momentum tensor first approaches the Allcn-Folacci value and only much later approaches the Bunch-Davies value. 
In Section VI we consider the geometric significance of the state independent asymptotic behavior of (T a b), relating 
it to the quantum effective action which determines the behavior of S c s un der global Weyl scaling and providing 
the generalization of the trace "anomaly" in the infrared, for £ ^ i . Section VII contains some discussion and final 
conclusions. There are two Appendices. Appendix A completes the proof of the Bunch-Davies attractor behavior in 
the cases of integer and pure imaginary u, while Appendix B contains a discussion of the simple harmonic oscillator in 
the limit of vanishing frequency, which shares many features with the spatially constant mode in de Sitter spacetime. 



II. SCALAR FIELD IN A RW BACKGROUND 



The metric for a general RW spacetime can be written in conformal time r\ in the form 



ds z =a z {r]) -dr] z + 



dr 2 



1 — KT' 



+ r 2 d& 



(2.1) 



Here a(i]) is the scale factor and k = 0, +1,— 1 corresponds to the cases of flat, spherical, and hyperbolic spatial 
sections, respectively. Throughout we use units such that h = c = 1 and the Misner, Thorne, and Wheeler pl| ] 
conventions for the curvature tensors, R% cd — T^ d c — ... and R a b — Racb- 

We consider in this paper a free quantum scalar field <I> with the quadratic action 



S = -\J d W^S [(V a $) 3 Qb (V b $)+m 2 <I> 2 +ei?<i> 2 ] 



(2.2) 



where V a denotes the covariant derivative, R is the scalar curvature, and g = dct(g a b). The mass m and curvature 
coupling £ are allowed to have any real value. The wave equation for $ obtained by varying this action is 



[- .09.09 



■eJZ]$(rj,x) = 



1 d 



a 4 dr] \ dr] 







— + m 2 + £B, 



$ = , 



(2.3) 



with A^ 3 ) the covariant spatial Laplacian. For spacetimes with the metric (2T) the field $ can be expanded as a 
mode sum in the form, M 



$(?7,x) 



a(rj) 



(2.4) 



where the integration measure is given by 



/d 3 k 



|d£(k) = j J^dkZ^ 



E 



k.l.m 



if K = , 
if K = -1 

if k = +1 



and the spatial part of the mode functions 5k (x) obeys the equation 

- A( 3 )y k (x) = (fc 2 - «)y k (x) , 



(2.5) 



with k = 1,2,... in the case of closed spatial sections, k = +1. The time dependent part of the mode functions tpk 
obeys the equation 



ipk" 



k' z + m 2 a 2 + I £ 



1 



a 2 R 



(2.6) 
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where primes denote derivatives with respect to the conformal time variable 77, and the scalar curvature in a general 
RW spacetime is given by 



R = 6 



(2.7) 



For the quantum field to satisfy the canonical commutation relations, the creation and annihilation operators are 
required to obey the commutation relations [ak, = ^ kk '' whereupon the ipk must obey the Wronskian condition 

Ml' ~ = i ■ (2.8) 
The components of the unrenormalized energy-momentum tensor (energy density and trace) are given by p2f 

1 



47r 2 a 4 



(6£-l) 



d M (fc)(2n fe + 1) {|Vfc| 2 + (fc 2 + m 2 a 2 )\^ k \ 



(2.9a) 



and 



£ u + 3p„ = (T) u = 



2Tr 2 a 4 



+ m - 1) 



d/j(fc)(2n fc + 1) <j -mV|^ fe | 2 + (6£ - 1) 

>2\ / 1 



-i^i 2 + -(^r+^) 

a 



k z + m'a z +[—-—]+[£-~]a'R 



(2.9b) 



where we have allowed an arbitrary number of particles in the initial state, n% = (a^ak), and the scalar measure d[i(k) 
is given by 



dfx(k) 




As we are considering spatially homogeneous and isotropic initial states (consistent with the RW symmetry), n k 
depends only on the magnitude k of the spatial wave vector k. Expectation values of the bilinears (ak<Zk) and (a k a ic) 
in a general state need not be considered since they can be removed by a time-independent Bogoliubov transformation 
at the initial time [Q] . Hence these initial state correlations may be parameterized instead by the initial data on the 
mode functions ipk ? together with the non- negative set of , with no loss of generality. 

Since the expectation value of the unrenormalized energy- momentum tensor (T a b)u is quartically divergent, a pro- 
cedure for defining a finite, renormalized expectation value must be given. We will follow the adiabatic regularization 
method p3|-p(| . In this method the renormalization counterterms are constructed using a fourth order WKB expan- 
sion for the mode functions. We denote these counterterms by (T a b)ad- They are given in Refs. (22| and p7| . The 
renormalized energy-momentum tensor is then 



(Tab^ren (Tab) u (Tab} ad 



(2.10) 



This subtraction scheme is not manifestly covariant in form, since space and time are treated quite differently. However, 
adiabatic regularization is equivalent to a covariant point splitting procedure in which the points are split only in the 
spacelike hypersurface of constant 77 |27],|28| , and the value of the renormalized (T a b) obtained by this procedure is the 
same as in a strictly covariant one. Hence this subtraction procedure does correspond to adjustment of counterterms 
to the quantum effective action, and (T a b)ren is covariantly conserved. As discussed in detail in Ref. p7| the adiabatic 
terms in all cases consist of an integral rather than a sum over k. The reason is that subtraction corresponds to purely 
local counterterms in the effective action, and thus must be independent of the global compactness or non-compactness 
of the spatial sections. 

For an arbitrary homogeneous and isotropic state to be p hysic ally admissible the renormalized energy-momentum 
tensor defined by the adiabatic order four subtractions in (2.10) must be both ultraviolet and infrared finite. In a 
general RW spacetime, ultraviolet finiteness for a fiel d w ith non-conformal coupling to the scalar curvature requires 
that the particular solution of the mode equation (2.6) in a general physical state must match the fourth order 
adiabatic form at large k, with the deviations from the fourth order WKB form falling faster than fc~ 4 . Likewise the 
initial number n k must fall faster than A;" 4 at large k, for the mode sums (integrals) to be ultraviolet convergent. This 
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is equivalent to the requirement that the two-point function of the scalar field have the vacuum Hadamard form |29| |32 
to sufficiently high order in the short distance expansion as the points app roach one another. As long as there are 
t wo linearly independent complex oscillatory solutions to the equation (|2.6|), the Wronskian normalization condition 



(2.8) can be imposed and the state will be free of any infrared di verg ences. However, when m 2 + £R — > for some 
low k in de Sitter space no such complex oscillatory solutions to ( |2.6| ) exist and the infrared finiteness requirement 
on the both the energy-momentum tensor and the two-point function of the physical state becomes non-trivial, as we 
discuss in detail in Section IV. 

A useful variation of the method of adiabatic regularization has been developed by two of us [j33f . In this method 
one first computes a quantity (T ao )d, obtained by expanding the adiabatic counterterms (T ao ) a d in inve rse powers of 
k and truncating at order k~ 3 . The same renormalized energy-momentum tensor defined in Eq. (2.10) is separated 



into the sum of two finite terms by adding and subtracting the simplified form of the divergent counterterms (T a b)d 

(Tab) ren — (Tab)n ~t~ (Talari 3 
(Tabjn = (Tabju — (Tab) d , 
(T a b)an = (Tab) d ~ (Tab) ad ■ (2-11) 

The full expressions for (T a b)d and (T a b)an are given in Ref. |3^] for a general RW spacetime. The advantage of this 
splitting is that (T a b)n and (T a b)an are separately conserved, and moreover, (T a b)an may be computed analytically in 
terms of the scale factor a(rf) and its derivatives ||33|| , Thus the state dependence of the renormalized (T a b) r en resides 
completely in (T a b) n , which can be computed numerically. 

III. THE BUNCH-DAVIES ATTRACTOR FOR < § 

We restrict our consideration henceforth to the particular maximally symmetric RW background of de Sitter space. 
The geometry of de Sitter spacetime can be described in a number of different coordinate systems. If k = the spatial 
sections are flat and the scale factor is 

(X 

a ( r i) = ■ ~°° < V < > k = , (3-1) 

V 

with a a real, positive constant, and R = 12a -2 . If n = +1 then the scale factor is 

a(rf) = a seer] , -~ < rj < ~ , K = +1 , (3.2) 

which is equivalent to a(ry) = a esc 77 with < rj < n. To simplify the notation we will generally use dimensionless 
units where a = 1 and R = 12, restoring the dimensionful quantities when it is instructive to do so. 

The asymptotic behavior of the energy-momentum tensor does not depend on k so the bulk of the analysis will be 
carried out in the flat (k = 0) coordinates. However, when we turn to the massless, minimally coupled limit in the 
next section, it will become useful to have a discrete k basis in order to separate out the k — 1 spatially homogeneous 
mode explicitly, since it is the most infrared sensitive. No confusion should be caused by our use of the same symbol 
rj for conformal time in both cases of flat and closed spatial sections, since we make use of only the k = coordinates 
in this Section and only the k = +1 coordinates in the next Section. We will not make use of the spatially open 
(k = — 1) coordinates i n th is paper. 

For the case of Eq. (|]l]) the general solution to the mode equation can be written as |l(]] [] 

MV) = \(-^e^ [dWH^i-kri+c^H^i-kij)] , (3.3) 

where the are Hankel functions and 

v 2 = ~ - m 2 a 2 - 12£ = - 7 2 . (3.4) 



x In Ref. (l6| the arguments of the Hankel functions are given as ki] rather than — kr\. We have chosen to use non-negative 
arguments to avoid complications that result from the fact that these functions have branch cuts along the negative real axis. 
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The latter notation is usefu l in the case v 2 < so that v = ij is purely imaginary. When v 2 > we will choose v to 



be the positive root of (|3.4| ) . From Eq. ( |3.3|) we see that solutions to the mode equation in de Sitter space depend on 
m and £ only through their dependence on the parameter v. Note that because of the minus sign in the arguments of 
the Hankel functions, it is the function Hu that corresponds to a positive frequency mode in the large k limit. The 



normalization of the mode function in (3.3) has been chosen so that the Wronskian condition ( |2.8| ) becomes simply 

(3.5) 

0) for all k. The renormalized 



|c 1 (*)| a -|c a (fe)| 3 = l. 
The Bunch-Davies (BD) state is defined by the choice, C\ = 1 and c 2 = (with n k 



energy-momentum tensor in the BD state is given by |15|,|1 



(Ta 



b) BD 



9ab 



64ir 2 



i? 



l , l 

R m 2 R - - 

18 2 



R z 



R 2 
2160 



- log 



/12m 



2 \ 1 



V R 



(3.6) 



where VK 2 ) = ^jgg£lg) j s the Higamma function. 

Th at this finite value of (T a b)BD coincides with the renormalized (T a fc) re „ defined by the adiabatic subtraction in 
( |2.10 ) follows from the fact that the BD state is an allowed fourth order adiabatic state. This may be checked by 
i 

comparing the asymptotic expansion of the exact BD mode function, ( -7r J') 2 {{^(—kri), for large k with the fourth 
order adiabatic mode function 

rn 



#(»?) 



1 



2W, 



(4) 



: exp 



(3.7) 



with W^'{rf) the fourth order adiabatic frequency. It is given explicitly in de Sitter space in flat conformal time 
coordinates by 



W^\v) = k + 



1 



1 



- v ~ - 



1 



1 



O 



2krj 2 \4 J 8fcV 

up to the required order at large k. 

For the general state with c 2 ^ to remain fourth order adiabatic, we must have for large values of k 

C(k) 



c 2 (k) = 



k 4 



(3.8) 



(3.9) 



for some complex function C(k) which vanishes in the limit k — > oo. This condition is necessary for an arbitrary 
(spatially h omog eneous) state to posses a finite energy-momentum tensor after the fourth order adiabatic subtraction 
defined by ( 2.1C ). Likewise the same condition of finite (T a b) requires us to restrict the average number of particles 
(a^ak) = n k by 

N(k) 



nk 



fc 4 



for some real function N(k) which vanishes in the limit k — > oo. The two ultraviolet conditions 

lim \C(k)\ = lim N(k) =0, 



(3.10) 



(3.H) 



on the physically allowed states guarantee that the Green's function for the scalar field is locally of the Hadamard 
form [p9|-p2[, and that the divergences of (T a b) match those of the fourth order adiabatic vacuum, and are re move d 
by the adiabatic subtraction procedure. We will call any state which satisfies the conditions ([j~9|), ( 3.10 ), and ( 3.11 ), 
together with the Wronskian condition ( |3.5| ), a UV admissible physical state ^|. 



2 In fact, the requirement of fourth order adiabatic states is slightly more restrictive than UV finiteness of the energy-momentum 
tensor in de Sitter space, since C(fc) can go to a non-zero cons tant at large k and the mode sums still converge. This is associated 
with the vanishing of the tensor, ^H a b, defined in Eq. ( |6.l[ ) in de Sitter space. If the field is non-conformally coupled, a state 
with C(k) —* constant would lead to a divergent energy-momentum tensor if the spacetime is not exactly de Sitter. If the field 
is conformally coupled and the spacetime is not exactly homogeneous and isotropic, then the energy-momentum tensor would 
again be divergent. Thus the most general physically acceptable UV states are fourth order adiabatic states. 
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We shall require also that the arbitrary physical state possess a two-point function and energy-momentum tensor 
which are free of any inf rare d diver genc es. Because the canonical dimension of (T a b) is four whereas that of (<$>(x)<f>(x')) 
is two, the conditions (3.9) and ( |3.1C| ) which require finiteness of (T a b) ren are more restrictive in the UV, whereas 



the condition of finiteness of ($(a;) < i>(x')) is more restrictive in the IR. These two sets of conditions will be sufficient 
to demonstrate that the energy-momentum tensor for any UV and IR admissible physical state approaches the BD 
value at late times for 3?(^) < |. 

To understand why such a result is to be expected and outline the mo re detailed proof which we give below, let us 
observe that at late times 7/ — > CP, the general state mode function (^^) behaves like 

~ (-»?)' """a" -4 - (3-12) 



Substituting this into (2.9a) and ( 2.9bj ) shows that to leading order at late times the contributions to the mode sums 



of the unrenormalized energy-momentum tensor behave like (— Ty) 3-2 " ~ a 2i/ ~ 3 for v real. Since the renormalization 
counterterms are state independent |2^| , the state dependent terms are the same in the unrenormalized and renor- 
malized energy-momentum tensor. One can perform all the UV renormalization in the BD state at a fixed time and 
collect the remaining finite state dependent terms which are unaffected by the subtraction procedure, and they all 
fall off at least as fast as (— J]) 3-2 " as rj — > 0~ for ^R(v) < |. 

These remaining finite state dependent terms in the energy density and trace are expressible as integrals over the 
wave number k with the general form 

I( V ) = ^ T R(k) S(-k V ) , (3.13) 

where R(k) is one of the four state dependent, but time independent functions 

|c 2 (fc)| 2 (l + 2n fc ), 

R[d(fc)c2(fc)](l + 2n fc ), 
9f[d(fc)c2(*)](l + 2n fc ), 

n k , (3.14) 

and S(z = —krj) is a product of the state independent Bessel functions and their derivatives. An explicit basis for 
the twelve products of Bessel functions Si{z) for i = 1, . . . , 12 which appear in th e integrals is given in Table 1. The 
essential point is that all the state dependent mode integrals of the form ( |3.13 ) are uniformly convergent for all rj 



(including r\ = 0) at both their lower limit, k = 0, and their upper limit k — oo, due to the IR and UV finiteness of 
the state. Hence the limit of r) — > 0~ can be taken inside the integral over k. Since, as Table 1 shows, all the Si(z) 
behave like 

S.t^s.o/', with A>0, (3.15) 
as z — > 0, for SR(i>) < §, it follows that 



r°° Ah 

lim J(r?) = / — R(k) Hm S(-kri) = 0, 

»7->0- Jo « »7-»0- 



(3.16) 



and all the state dependent contributions to (T a fc) ren vanish at late times. 

The validity of bringing the limit i nside the integral depends on the uniform convergence of the integral at both its 



upper and lower limits. In the form (3.13) the behavior of the S(—krf) factor at small arguments (and the absence of 



any IR divergence from the R(k) factor) clearly guara ntees the uniform convergence at the lower limit. However, the 



change of variables z — ~kr\ brings the integral ( 3.13|) into the form 



m = r t r (^) s(z) - (3 - i7) 

In this form it is clear that the uniform convergence of the integral at its upper limit is guaranteed by the falloff of 
the state dependent mode functions at large k, namely 

lim R(k) = lim R ( — ) = . (3.18) 
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Equations ( [3.15 ) and (3.18) guarantee that both the IR and UV contributions go to zero as r\ — > for 3t(u) < |. 
In fact, we can go one step further by using the information from the fourth order adiabatic nature of the state 



lim k A R(k) = , 



(3.19) 



to conclude that the UV contribution to the integral from z ~ 1, and very large k ~ (— rj)^ 1 in (3.17) falls faster than 
(— r)) 4 as rj — > 0~. This is faster than the IR contribution which falls only as {—rj)"* if /3j < 4. If /3, > 4 then both 
the IR and UV contributions fall faster than (— ?y) 4 . Hence, the Si(z) for which /3j > 4 are subdominant at late times. 
Thus we conclude that the leading order state dependent terms at late times are those with the smallest These 
give an IR dominant, (i.e. finite t, z « 1) contribution of the form, 



/>oo 

7(77) -» (-'?) ,3! / d/fc A^" 1 i?(fc) -» (- 
Jo 



dk k 



2-2v 



R(k), 



(3.20) 



as rj — ► . This last integral is guaranteed to converge at its upper limit for all 5ft(^) > and in particular for 
< %t(v) < | by (3.19). Hence for all 9?(f) < | the state dependent contributions go to zero as a 2 "" 3 at late times. 



The limiting case when v 



is special because then the state dependent IR contributions apparently does not go 



to zero at late times. This case will be considered separately in the next Section. 

We will now make the proof more explicit by giving the form of all of the state dependent ter ms o f the scalar field 
in de Sitter space, and analyzing the IR and UV contributions in detail. If we make use of Eq. (3.5), we find that in 
an arbitrary state 



{Tabjren — (Tab) BD + (Tab) 



b/SD 



(3.21) 



where (T ao )sD is the finite state dependent term, depending on the coefficients ci(fc), C2(k), and n k , which may be 
expressed as an integral over the wave number k in the form 



(Ta 



b/SD 



1 

4^2 



dk Iab(k,T)) . 



(3.22) 



The explicit expressions for the integrand I a b depend on whether v is real and not an integer, real and an integer, 
or pure imaginary, although the result will be the same for all < |. In the rest of this Section we restrict our 
discussion to the case where v is real and not equal to an integer. The cases of integer and imaginary values of v are 
covered in Appendix A. For real v after some regrouping of terms we have 



I°o = A 1 (k) 
+A 2 (k) 
+A 3 (k) 
/ = Ai(A) 
+A 2 (k) 
+A 3 (k) 



Sx + i(9- 
S2 + \(9 



Am 2 



48£)S 4 --(4£-l)S 7 + Si 



S 3 + -(9 



Am 



Am 



•480S 6 --(4£- 



48£)S 6 -~(4£. 



-2(6£-l)Si--(9 + 8(3£-l)m 2 



1)5*8 + Sn 
l)iSg + Si 2 , and 
- 102£ + 288£ 2 ) Si + 3(6£ - 1)S 7 



(3.23a) 



-2(6£-l)S 2 --(9- 
-2(6£-l)S 3 -i(9- 



8(3£ - l)m 2 - 102£ + 288£ 2 ) S 5 



■ 8(3£ - 1)?ti 2 - 102£ + 288£ 2 ) S 6 



2(6£-l)Si 
3(6£-l)S 8 + 2(6£-l)Sn 
3(6£- l)S 9 + 2(6£- l)Si2 



(3.23b) 



with 



CSC 



'(vtt) ((1 + 2n k )\c 2 \ 2 + n k ) + - (l - cot 2 (^7r)) (1 + 2n k ){ Cl c* + c*c 2 ) 



+— cot(z/7r)(l + 2n k )(c\c* 2 — c*c 2 ) 



(3.24a) 



A 2 (k) 



cot(^7r) csc(^7r) ((1 + 2n k )\c 2 \ 2 + n k ) cot(i/7r) csc(^7r)(l + 2n k )(ciC2 + c{oi) 
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+ - csc(^7r)(l + 2n k ){c\c* 2 - c^c 2 ) 



csc 2 (^tt) ((1 + 2n k )\c 2 \ 2 +n k ) - -csc 2 (^tt)(1 + 2n k )(c lC * 2 + c\c 2 ) 



(3.24b) 
(3.24c) 



and the Si{—krj) composed of various products of J v {—krj), J^ u (—kr)) and their derivatives. Making use of the general 
formula for the product of two Bessel functions pm , 



p=0 



(-lF(§) W2p I> + / , + 2 P +l) 



(3.25) 



we can expand the Si in power series of the form, 



Si (z) — S ' Si t p 



(3.26) 



p=0 



with 



-krj. The explicit expressions for [3i and S% for the case of real, non-integer v are given in Table 1. 



1 


Qi 

HI 


Si(z) 


1 


h + 2v 


z^Jl(z) 


2 


5 


z 5 J l/ (z)J- l/ (z) 


3 


b-2v 




4 


3 + 2v 




5 


3 


z z J u (z)J- u (z) 


6 


3-2^ 




7 


3 + 2^ 




8 


3 


z 4 A (J V ( Z )J_ U ( Z )) 


9 


3-2^ 




10 


3 + 2^ 




11 


3 




12 


3-2v 





Table 1 



We are interested in the behavior of the finite state depend ent te rms, (T a b)sD, in the limit 77 — ^ .To investigate 
this limit in detail it is useful to break up the integral in Eq. (3.22) into the three parts 



(T a b)sD = ^j— j / dfc I a b(k, 77) + — — / dfc I a b(k, rj) + j dk I a 



.{k,rj) 



(3.27) 



Here A is a finite positive constant. For k > A we can make use of the UV conditions (3J5), ( 3.1C ), and ( 3.11 ) in 
the second and third integrals. Thus, the most infrared sensitive integral is the first one. The positive constant Z is 
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arbitrary, provided only that Z > —Xt], which is always satisfied for fixed A and small enough — r\. Hence the second 
integral provides the bulk of the contribution of the state dependent wave numbers that have redshifted outside the 
de Sitter horizon at late times. The third integral is the contribution of the state dependent terms still within the 
horizon which go to zero very rapidly at late times due to the UV conditions. If the 77 — ► 0~ limit is uniform then all 
three integrals should vanish unconditionally in this limit, i.e. without any restrictions on the arbitrary parameters 
A and Z. 

We begin by analyzing the first integral in (3.27). We are considering only fourth order adiabatic states that are 
IR admissible. Hence the fc integration converges at its lower limit and we can expand the integrand for this first 



integral in a series of the form ( |3.26 ) and interchange the order of summation and integration. Each term in the 
resulting sums contains an integral over fc that is finite by assumption and a factor of (— ^) 2 P+ft where p = 0, 1, 2 . . . 
is a non-negative integer. In fact, since at late times, — A?y <C 1 for any finite A, it is sufficient to consider only the 
leading p = term. From Table 1 it is clear that for v < |, > for all i. Hence in the limit 77 — > — , the integrand 
vanishes and therefore, the fi rst in tegral on the right hand side of Eq. ( |3.27 ) vanishes in this late time limit. 

For the second integral in (3.27) we may utilize the expansion (3.25) again. Since the integration is between finite 
limits for finite 77 we can exchange the order of summation over p and integration over k. Having done so we can then 
bound each term in the sums by taking the absolute value of its factors. Inspection of the expressions (3.23t) and 
( |3. 24c ) indicates that the result is a linear combination of terms involving integrals of the three possible forms 



Zi = 
Ti = 
1 3 = 



z ^dk 
k 

- z/r > dfc 
- z/r > dfc 

T 



|c 2 (fc)| 2 (1 + 2n k ) (-fcry) 2 ^ , 
|ci(fc)c 2 (fc)| (1 + 2n k ) (-k V ) 2 ^ 
n k (-kr,) 2 ^ , 



multiplied by constant coefficients. Because of (3.11) for any A 



|C(fc)| < C, 
N(k) < N , 



(3.28) 



(3.29) 



for some real pos itiv e numbers C and N, depending on A. With these bounds we can bound the values of |ca(fc) j and 
n k , and use Eq. (3J5) to bound |ci(fc)| for all fc > A as follows 

C 



Mfc)| < 



fc 4 



c 2 



c 2 



1 < |ci(fc)| = [1 + |c 2 (fc)| 2 ] 3 < 1 + |c 2 (fc)| 2 < 1 + < 1 + — , 



A 8 



n k < 



N 

fc 1 



2N 

l + 2n k <l + j r <l 



2N 
A^ 



(3.30) 



Using these bounds the integrals in Eq. ( 3.28 ) can be bounded as follows 

-z/r, 



Tj < C 2 (l + 2iV/A 4 ) J dfc k 2p+l3 '- 9 (-T 1 ) 2p+i:! ' 



C 2 (1 + 2N/X 4 ) 
fa + 2 P - 8 



[(- V ) 8 Z 2 P + ^ 8 - (_ 7? )2p+/3, A 2 P +ft~8] 



1 2 < C(l + C 2 /A 8 )(l + 27V/ A 4 ) 



-Z/r, 



dfcfc 2p+ft - 5 (-r/) 2 f+ ft 



C(l + C 2 /A 8 )(l + 2N/X i ) r/_ s4^2p+/3i-4 _ r \2p+ft iiSp+ft 

ft + 2p-4 LV " K " 



J 3 < N J Zh dkk 2 P+^ 5 (-r 1 ) 2 P+^ 
N 



ft + 2p-4 



(3.31) 
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Each of these bounds vanishes in the limit 77 — ^ . Therefore each of the terms appearing in the second integral in 
( 3.27 ) vanishes in the late time limit for 5f(i/) < §. 

If v = I some of the terms will have vanishing denominators and should be interpreted according to the limiting 
relation 



lim 



(-1)' 



- X q 



Q 



= log 



-77A 



(3.32) 



but the appearance of these logarithms does not change the result since they are always multiplied by at least 
(„ 7/ )2 P +ft which van i s hes for 5ft(z/) < | for p = 0, 1, 2 . . .. The cases v = 0, 1 also involve logarithms in the Bessel 
function expansions but the result that the second integral in ( |3.27 ) vanishes in the late time limit is unchanged. One 
may consider also the case when v — 77 is pure imaginary, where the forms of the coefficients (3.24c) and bilinears 
Si(z) in the table change somewhat, with again the same result. For completeness these cases are treated in detail in 
Appendix A. 



Finally, for the third integral on the right hand side of Eq. ( 3.27 ) we do not expand the Bessel functions in powers 
of (— krj). Instead we change the integration variable to z = —krj and find integrals of the forms 



Ji = (-vf 



dz 



C 



-T) 



N 



J 2 = (-rj f I _ (3? or 3) <j Cl ( — 

Js = (-vTI 5«[-|S,( Z 



C* I — 

V 



Si(z) 



1 + 21-2} N IJL 



Si(z) . 



(3.33) 



Since we are considering fourth order adiabatic states these integrals are finite for all —00 < rj < 0. Furthermore, the 
integrands of these integrals are finite throughout the entire integration range for — 00 < 77 < 0. With the lower limit 
Z > fixed, it is clear that we can evaluate the integrals at 77 = by first taking the limit 77 — > 0~ of the integrands 
and then computing t he in tegrals. Since the integrands all vanish in the limit 77 — > — , the integrals do as well. Then 
the third integral in (3.27) vanishes in the late time limit for 3?(^) < |. From the adiabatic four conditions (3.11) 
these UV contributions go to zero faster than (— ?y) 4 for any fourth o rder adiabatic state. 

Therefore we have proven that for 3?(^) < | all the integrals in Eq. ( |3.27 ) vanish in the limit 77 — * 0" for an arbitrary 
state that is both infrared and ultraviolet finite. The fact that all three integrals vanish at late times, independently 
of the parameters A and Z which we introduced to control the IR and UV contributions of the mode integral verifies 
that the total mode integral does converge uniformly in 77, as expected. Hence all the state dependent contributions 
vanish at late times, 77 — > - , and we have proven that in this limit (T a b) ren — > (T a fc)siJ- Thus, we conclude that 
for tft(v) < I the value of the energy-momentum tensor in any physically admissable homogeneous and isotropic RW 
quantum state asymptotically a ppro aches the Bunch-Davies value in de Sitter space at late times. 

Moreover, from inspection of ( p.3l| ) we observe that the contribution from the upper limit at Z in the second integral 
always falls faster than (— 7i) 4 , while the contribution from the lower limit at A falls off only as (—77)^ as 77 — * CP. 
This must be the same order as the first IR integral so that the arbitrary parameter A drops out of the final result. 
Hence our detailed evaluation has verified that the leading order state dependent corrections to the BD expectation 
value at late times come from the terms in (T , a b) ren with the smallest /3,-. From Table 1 these are the i = 6,9, and 
12 terms. Collec ting t hese terms from ( 3.23b[ ) and (3.24c), and the numerical coefficients from the Bessel function 
product formula (3.25), we conclude that the leading order behavior of the energy density and trace at late times are 
given by 



e — > a 



T — » 2a 



2i/-3 



!- in- - -'Hi; ■ 1 !(/' ■ 1 1 



poo 

/ dfc k 2 - 2u R{k) , and 
Jo 



-m 2 + (6£ - 1) 



I'--) +2^ + m^ + 2(6^-l) 



l jf dfc k 2 - 2v R{k) . 



respectively, with 



R{k) 



2 2t/ - 3 csc 2 (7r^) 
k[T{-v + 1)] 2 



{(l + 277 fc ) [\c 2 (k)\ 2 -n( Cl c* 2 )] +n k } 



(3.34a) 
(3.34b) 

(3.35) 
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The leading order state dependent contribution at late times is a 2v 3 from the IR part of the mode integral, which 
falls off very slowly if Ji(z/) — * |. We examine this latter limit in detail in the next Section. 



IV. THE ALLEN-FOLACCI ATTRACTOR FOR 



In the analysis of the previous Section we saw that the terms with the slowest falloff at late times were those 
with the smallest (3i — 3 — 2v which behave like a 2lJ ~ 3 for v real and positive, and the coefficient of this falloff is 
controlled by the finite k, IR part of the mode integral. To examine the limit v — > § carefully, it is easiest to work 
with closed spatial sections and a discrete set of mode functions in order to treat the most infrared sensitive, spatially 
homogeneous k = 1 mode separately from the rest, instead of dealing with an infrared sensitive continuous mode 
integral. 



T he scale factor for k = +1 is given by Eq. (3.2) and under the variable substitution Q — i tan 77 the mode equation 
( [2.6|) becomes Legendre's differential equation. Hence the general solutions may be expressed in terms of associated 
Legendre functions P^i (C)- Since as conventionally defined these functions have a cut discontinuity on the real 
axis from — 1 to 1 if k is an odd integer, we write the fundamental complex valued solution for real v in the form 



fk(v) = 




(i tan 77) 



v;k + l 



1 — i tan 77 



(4.1) 



where e(rj) = d(rj) — d{—rf) = ±1 is the sign function and F is the hypergeometric function. The phase factor depending 



on the sign of 77 for odd k removes the discontinuity in the P_i, function as 77 approaches zero from positive or 



negative values, as the second form of (4.1) makes clear, since the F function is an analytic function of ^rrS with no 



branch cuts for £ on the imaginary axis. With the normalization factors chosen as in (4.1) the Wronskian condition 



fkfk ~~ fkf'k — ' ■ 



is satisfied. Hence the general solution of the mode equation (2.6) is 
with 



(4.2) 



(4.3) 



l/3 fe | 2 = 1 



(4.4) 



The Bunch-Davies state is given by = 1 and (3^ = 0. 

Now as v — ► I inspection of ( f4.1[ ) shows that all the fk for fc > 1 are regular. In fact, in that limit the hypergeometric 
series for F terminates and the fk{r/) become the elementary functions 



lim f k (t]) 



-ikr/ 



[2k(k 2 - 1)]- 



(k + i tan 77) . 



2,3,. 



(4.5) 



so they can be treated in the Bunch-Davies state au = l,/3fe = with no difficulty. However in this limit the k = 1 
mode function is singular and must be treated separately. The behavior of the k = 1 mode as v — > | is similar to 
that of a simple harmonic oscillator mode as its frequency goes to zero, i.e. in the limit where the harmonic oscillator 
becomes a free particle. The zero frequency limit in this simple flat space analogy is reviewed in Appendix B. Just 
as in that case, one can construct regular solutions to the mode equation in the limit v — > | by taking suitable linear 
combinations of f\ and /*. In fact, the limiting form of can be found from Sec. 2.3.1 of Ref. p5|, which gives 



sec 77 



- yi - v scc v (v 



sin 77 cos 77) 



(4.6) 



We have neglected terms in the real part of f\ that are of order J | — v. We have also neglected all terms that go to 
zero faster than this as v — 



Extracting the scale factor 0(77) = sec 77 we now define the real functions u and v by 
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u(rj) 



We next define new coefficients 




(4.7) 



A 



-i lim 



- v {ax - Pi) 



B = lim 



2,/| 



which arc finite in this limit. With these definitions the normalization is such that 

A*B — B*A = i , 

and the limit of the general k — 1 mode function may be written 

'A 



lim ipi(r)) = sec 77 (Av + Bu) = a{rf) 

One can also define the time-independent Hermitian operators 

1 



(77 + sin 77 cos rj) + B 



(at +/3i) 01 + (ai +/3i)* a| -> Pai + P*aJ 



P = — m/ 

2 



{a x - 01 - (ai - /3i)* a| -> A01 + A*aj , 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



obeying the canonical commutation relations, [Q, P] — i. They also remain finite in the limit v — > |. 

We will consider the limit v — > | of the energy density in the particular order of fixing the mass of the scalar field to 
be m = and letting the dimcnsionless parameter £ approach zero, since this case is relevant for the infrared scaling 
analysis of Section VI. From (3.4) with m — and £ small we have 



(3 - 2u) 
8 



(4.12) 



The complementary case, £ = and m 2 — > is similar and has been discussed in Refs. [[[7) and p9[ . From ( [2.9a ) we 
read the energy density in the k = 1 mode for 771 = 

si = -<T° )*=i = {Kl 2 + l^i I 2 + (<*- l)[tanTj + V>iVi) " (tan 2 77 - l)^ 2 ]} . (4.13) 

We are interested first in the asymptotic form of this e nergy density at late times, 77 — > J, and then in the limiting 
form of the resulting expression as £ — > according to ( 4.12 ). The asymptotic late time limit of the mode function 
/1 (77) for any v > can be found from the inversion transformation of the hypergeometric function, given by formula 
2.1.4 (17) of §|. We find that 



2r 



/I 



r(| + ^) 



H) 



i sec 77 



sec 77 /1 ~ a y+ 2 



(4.14) 
(4.15) 



as 77 — > which is the same late time behavior in terms of the scale factor that we found in the spatially flat 
coordinates. Notice that at v = § the phase factor cancels and these limiting forms of the oscillatory mode functions 
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become real, while the T function has a pole singularity there. Thus, keeping only the leading behavior as v approaches 
| we find that the late time limits of the mode function and its derivative are 



ai + (3\ 



(4.16) 



Thus, the dominant terms in the k = 1 energy density e\ in this limit are |V'1| 2 and the terms involving = tan 77, 
and we obtain 



lim E\ 



3(1 + 2n x ) 



32tt 2 



|ai+/3i| 2 + C(3-2z/) 



,2i/-3 



(4.17) 



The singularity in the T function has canceled against the £ in the numerator. The remaining coefficient is clearly 
state dependent. From this asymptotic form of the energy density in the k = 1 mode at late times it is clear that for 
v close to but still less than |, the state-dependent energy density E\ goes to zero at late times, albeit very slowly, 
which is consistent with our previous flat section analysis. 

The corresponding expression for the trace in the k = 1 mode is 

- ei + 3 Pl = (1 2 ^ 2 2 Q 4 l} (6e - 1) {-Kl 2 + [tan r? (^f + + [sec 2 r, + 2(6£ - 1) sec 2 V Mi\ 2 }} ■ 



Substituting the late time asymptotic forms (4.16) to leading order in £ as before yields 

3(1 + 2m). 



ei + 3pi 



-|ai+/3i| 2 + 0(3-2i/) 



Jlv-Z 



-4ei 



(4.18) 



(4.19) 



Hence, p\ — > —ex and the contribution from this mode is de Sitter invariant at late times for any initial physical state. 

Since the Bunch-Davies state is au = 1 and j3u = — for all k an d this state has a renormalized (T a b) which is 
strictly time independent for all v, the time dependent contribution in (4.17) and (4.19) for v < | in the k — 1 mode 
must be canceled by a time dependent contribution from all the other modes in the renormalized energy-momentum 
tensor. In other words, the late time behavior of all the k > 1 modes in the BD state must be 



£bd 



k>l 



£bd — 



327T 



+ 0(3-2u) 



(4.20) 



for v close to § . The pressure for the k > 1 modes is obtained from this by the de Sitter invariant relation p 



The subtraction of the second term in (4.20) can be understood in a different way. Consider the short distance 



expansion of the BD two-point function for v < | , namely |L9l 



G 



L]D (x,x')^^ 



— - log(l - Z) + 3^- 



1 - z 



(4.21) 



for the de Sitter invariant bi-scalar, Z(x,x') — > 1 and v close to |. The constant term is singular at v = |, but 
it gives a finite contribution to the energy density from the £,G a b4> 2 term in the energy- momentum tensor. In the 
dimensionless units we are using this is equal to 



— i (- 

8tt 2 ? I 12 



1 



32tt 2 



(4.22) 



for £ — > according to ( 4.12 ). The last constant term in ( 4.21 ) is absent in the short distance expansion of the 
Allen- Folacci two-point function JD]], and in the Bunch-Davies two-point function it comes entirely from the k = 1 
mode in the BD state. Hence the contribution to the e nergy density of the k > 1 modes in the BD stat e does not 
contain (4.22), and must equal the full Ebd minus (4.22) from the k = 1 mode, which is equivalent to (4.20) at v = |. 

With the k = 1 and k > 1 mode contributions separated we ar e now in a position to take the limi t of v — > |. 
Apparently we would obtain a state dependent contribution from ( 4.17 ). However, inspection of fl4.10| ) shows that 



14 



requiring the k — 1 m ode function to be regular as v 
this limit. From (4J3) we then have 



is equivalent to requiring that A and B remain finite in 



ai + 01 



2B 



0. 



(4.23) 



Therefore the apparently state dependent contribution (4.17) from the k = 1 mode van ishes in any infrared finite 
state parameterized by finite A and B, while the k > 1 contribution is still given by ( 4.20 ), and we have the result 



e R = lim Ebd (0 



m=0 



32tt 2 



777 = ^AF 



(4.24) 



in any (A, £?) "vacuum" state with oik = 1> /3& = nk = for fc > 1, upon restoring physical units. The difference here 
is just that required to give the Allen-Folacci (AF) renormalized exp ectat ion value, pro vided the limit of the BD value 
is taken in the same order of m = 0, £ — * that we have evaluated ( 1.17| ). Since from (3.6) 



£bd(Q 



m=o I6n 2 

and both the BD and AF values are de Sitter invariant 



180 6 V s ' 



(4.25) 



lim (T ab ) ren (£) 



m=0 



3gab 
167T 2 

119 

960vr 2 



R 
12 

R 
12 



9ab 



1 1 1 

180 ~ 6 ~ 2 
119 R 2 



138 240tt ; 



■9ab — {Tab) 



b AF 



(4.26) 



and we obtain the de Sitter invariant AF result for all (A, B) "vacuum" states in the massless, minimally coupled 
case. 

With the above careful analysis of the spatially homogeneous mode and the independence of the asymptotic value 
of its energy- momentum tensor in an arbitrary k = 1 infrared finite state characterized by A and B, it is now 
straightforward to carry out the proof of the attractor behavior of the AF state for an arbitrary UV fini te p hysical 
state with m = £ = 0, by all o wing t he k > 1 modes to have 0k and n k different from zero. Substituting ( [4. 3D , (fOj), 
( |4.10|) , and ( [4.161) into fl2.9a|) , (|2.9b[ ) with m = £ = 0, and using (Q gives 



■{T° ) AF + {l + 2n 1 ) 



|^4| 2 cos 6 ?7 



1 °° f 1 
4-^2 X, 1 t 2nfc + + 2n fe)l/ ? fc| 2 ] [ fc3 cos4 V + k (~ cos4 V + 2 cos2 »?)] 



(l + 2n fc ) 



(Aa;e 2,fc " + 0* k a k e- 2lk, i)k{- cos 4 + 1 cos 2 n) 



+ i (0 k a* k e 2ikr > - 0* k a k e^ 2lkr ')k 2 cos 3 ry sin 77] } 
{T) ren = {T)af + (1 + 2m) ^ 

4 -^£{[2n fc + 2(l + 2n fe )|/3 fc | 2 ] fccos 2 r; 



(4.27a) 



k=2 



(1 + 2n k ) [{0 k a* k e 2lkr < + 0* k a k e- 2l ^)(-2k 3 cos 4 r) + k cos 2 77) 
- 27 (0 k a* k e 2 ^ - P* k a k e- 2 ^)k 



2 cos 3 n sin 77] } 



(4.27b) 



cos 2 77, 



Provided the k sums converge, it is clear that all the state dependent terms contain at least one factor of a" 
and so vanish in the limit of 77 — > However, the requirement that the state be fourth order adiabatic just guarantees 
this convergence, for the same reason as in the previous analysis in spatially flat coordinates. Indeed we have 



\0k 



n k 



C{k) 

fc 4 ■ 
N(k) 
~k*~ 



(4.28) 
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for some C(k) and N(k) that vanish as k — > oo. This is sufficient to guarantee the absolute convergence of all terms 
in the sums. Since all state dependent terms are multiplied by at least two powers of cos 77 = a -1 , which vanishes 
in the late time limit rj — > we conclude that any UV and IR admissible state of the massless, minimally coupled 
scalar field has an energy-momentum tensor which approaches the AF value, (T a b)AF in the late time limit 77 — > 

If one considers the contributions of the state dependent terms to the energy density and trace from the k > 1 
modes, for v not exactly §, the kinematics is essentially the same as our previous analysis of the k — 1 mode. 
Their contribution also falls off proportional to £,\ak + /3fc| 2 a 2l/_3 at late times for v close to §. However, there is no 
compensating large factor coming from the pole in the T function normalization constant as there is for the k — 1 
mode. Hence, the coefficient of this slow fall off goes to zero as the massless, minimally coupled limit is approached. 
That is, exactly at m = £ = when the contribution of the k > 1 modes no longer falls off and they could in principle 
contribute to the asymptotic value of {T a b) re n at late times, at that very point their coefficient vanishes identically 
and they make no contribution at all. Thus, at precisely m — and £ = the difference between the BD and AF 
values can be attributed entirely to the additional condensate in the spatially homogeneous k = 1 mode alone, and 
there are no slow transient modes in our explicit analysis of the massless, minimally coupled energy and trace. 

These considerations are relevant to the case when v = | but m and £ are separately non-zero. The calculation 
is almost identical but the conclusion is different, since now the finite k > 1 mode sum is multiplied by a coefficient 
£ which does not vanish. The entire mode sum of state independent terms from k = 2 up to a{rj) do not fall off at 
late times and in fact all add up, to give a contribution proportional to ($ 2 ) ~ X)fc=2 ^ ~ l°g a m (T g i,) rew , which 
grows linearly in comoving time. The explicit expression is most conveniently calculated by using (2.11) to divide 
the energy- momentum tensor into a state dependent numerical part, (T a b) n and a state independent analytic part, 
{T a b) a n- T hey ar e se parately co nserv ed 1 33 . Th e quantities {T ao )d and (T a b) a n are given in Ref. p3[ . By substituting 
Eqs. ( ff .16| ) and ( El] ) into Eqs. (2.9a) and ( |2.9b| ), and subtracting the expression for (T a b)d given in Ref. (3^], we find 
that {T a b) n approaches a state dependent, finite constant in the limit ry — > 5. However, the quantity (T a b)an has a 
term that is proportional to (3> 2 ). At late times this term behaves as log a and dominates. In fact 



3£ 

£ren ~> - ^~ J log(cOS 77) + £ft , 

(T) r en ~> "o log(cOS f]) + q 2 , 
TX 1 



(4.29a) 
(4.29b) 



with the constants q\ and qi dependent on the state of the field and constrained by the conservation equation. Since 
log(cos 77) — >• i in comoving time, the energy-momentum tensor grows linearly with t at late times. For £ > (and 
to 2 < 0) the linear growth in comoving time decreases the effective cosmological "constant" , whereas for £ < (and 
to 2 > 0) it increases it. In either case the back-reaction of the energy-momentum of the quantum scalar field for 
m 2 + £i? = (but m 2 = — £i? ^ 0) on the geometry certainly cannot be neglected at late times since (T a b) grows 
without bound for any physical state. The fact that a massive non-minimal field with £ < can induce an effective 
cosmological "constant" due to inflationary particle production was noted in Ref. using a different approach. It 
is an interesting open question whether this linearly growing behavior (in proper time) carries over to the physically 
more relevant case of one-loop gravitons, since the mode functions for gravitons in a particular gauge obey the same 
equation in a RW spacetime as do the mode functions for a massless minimally coupled scalar field |37|. 

Finally, we note that for all real values of v, the analysis of the k = 1 mode in this Section may be extended to all 
of the highe r k m odes, since the late time behavior of the higher k modes is determined by that of the hypergeometric 
function in (4.1), which gives 



r(fc 



2r k- 



T(2u) (~i) k ( isecr] 



(4.30) 



Thus all of the higher k modes behave as a v 2 and give the (unrenormalized) energy density the leading order late 
time behavior 



v 

2 



12^-1) 



F(2^) 



2T (*+!) 



] 2 f> + 2 "*) r ^j^l> * + e**"""-^ ■ (4-31) 



v) (fc!) 



which grows (or shrinks) like a 2v 3 for an arbitrary physical state unless either the coefficient vanishes or + 
e i7r(fe+!/-i-) _ q £ Qr a j| y. However, the latter is impossible since it is inconsistent with the requirement that all 
physically allowable states be fourth order adiabatic states, which requires that ak — > 1 and (3k — * sufficiently fast 
as k — > 00, in order for the renormalized energy density to be UV finite. Since the fourth order adiabatic subtraction 



1G 



behaves at most like a constant at late times 
the sum in ( 4.31 ) converges at large k 



the a 2v 3 behavior is not affected by the UV subtraction, and indeed 
for any UV admissable state. Hence after the adiabatic four subtraction the 



leading order late time behavior indicated by (4-31) survives unles s the f irst bracket in front of the entire expression 
vanishes. The quantity in this bracket is identical to the factor in ( 3.34a ) found in the flat spatial secti on an alysis of 
the last Section. The corresponding quantity for the trace is given by the factor in curly brackets in (3.34b), which 
is a similar combination of v, £, and m 2 . If we require that both of these factors vanish identically, to eli mina te the 
leading order behavior in all components of (T a b}ren, then these two conditions plus the defining relation (|3.4|) give 



in 



v(2v-Z){2v- 1) 
4(^-2) 
Q - 3) 
8(1/ - 2) ' 



(4.32) 



Thus, except for v = 2, for any given v there is always one value of m 2 and one value of £ for which the coefficients of 
these leading order terms in (T a b) ren vanish. The next to leading order terms go like o 2l/_5 , and grow without bound 
at late times in any case when v > |. 

The above analysis implies that, for most values of m 2 and £, when v > | the leading order terms in the components 
of the energy-momentum tensor grow without bound like a 2u ~ 3 in de Sitter space for any physically admissable initial 
state of the scalar field. These values of v correspond to the purely tachyonic cases m 2 + < 0. 



V. NUMERICAL STUDIES 



In this Section we display numerical results for various values of v. All results are given in dimensionless units 
where a = 1 and R — 12. The quantum state in each case is a fourth order adiabatic state matched to the vacuum 
at some initial time 770. That is, we choose initial conditions for the mode function and its first derivative to be 

-j— Wo) = — f-(Vo), (5-1) 



with ip^fi the fourth order adiabatic mode function defined by (3.7) with phase measured from the initial time tjq. 



The initial time was chosen to be to = 1 in comoving coordinates, i.e. secijo — cosh(l) = 1.54308 . 



For 3?(f) < I the proof in Section [II states that the energy-momentum tensor approaches the Bunch-Davies value 



(3.6) at late times for an arbitrary physically admissable state. This occurs due to a redshifting of the state dependent 



part of the energy-momentum tensor. The initial state dependent transient contributions fall off like 

a 2 — 3 = (cosht) 21 '- 3 . (5.2) 
Thus the characteristic time to approach the BD value is 

' (5.3) 



3-2^ 



In Fig. 1 we plot the renormalized energy density for these adiabatic initial conditions with m = and £ = y. For 
this relatively large value of £, the characteristic time r is of order one and the energy density approaches the BD 
value within one expansion time. 

For smaller values of £ the initial value transients persist for longer times. Our analysis in the previous Section 
shows that the k = 1 contribution is essential for the shift from the BD to AF value as £ — > 0. Since the k = 1 mode 



in an arbitrary physical state contributes to the energy density the value £1 given by ( 4.17 ), up to terms which fall 
off like a~ 2 , for arbitrary ai,0\, and m, while the BD state has otx = = n\ = 0, it is clear that the difference of 
the renormalized energy density from the BD value coming from this mode is 



R- 2 



For the adiabatic initial conditions here we find that 



17 



(ai +/3i)sec?7o 



(5.5) 



Since W^M remains finite as £ — > 0, provided that 770 7^ 0, it follows that in this limit \a\ + /3i| 2 ~ 3 — 2f <~ £ also 
goes to zero. Hence for small £, goes to a value close to ( but slightly larger than) the AF value after a time of 
order one. This is observed in both Figs. 2 and 3 for £ = and £ = j^q, respectively. 
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FIG. 1. The energy density for m = and £ = i as a function of comoving time t for adiabatic initial conditions at t = 1. 
The attractor behavior of the Bunch-Davies value for the energy density is illustrated. 
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FIG. 2. The energy density for m = and £ = yij as a function of comoving time t for adiabatic initial conditions at t = 1. 
The line segments at right giving the values of the energy density for the Allen-Folacci (AF) state, and the Bunch-Davies (BD) 
state (for m = and £ = ) . 
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FIG. 3. The energy density for m = and £ = for the same initial conditions as the previous figures. Notice the very 
long-lived transient behavior of e, staying close to the AF value with a small positive linear slope towards the BD value. 



Further, the energy-momentum tensor contains a term proportional to ^($ 2 ) which would grow linearly in comoving 
time for v very close to I, except for the factor of a 2l/ ~ 3 that damps it to zero at very late times t>r, Hence on 
times 1 < i « r, where the a 2ly_3 factor is essentially constant, we should expect this linear growth of e in comoving 



time with a slope proportional to £, given by (4.29a). This behavior is demonstrated in Figs. 2 and 3. In the latter 
case £ is so small that a\ + (3\ and the slope are nearly vanishing and the energy density stays close to the AF value 
until times of order r = 125, which is much larger than the times shown in Fig. 3. When £ — (still keeping m = 0) 
both ol\ + [3\ and the slope vanishes identically, so the energy-momentum tensor goes to the AF value and remains 
there. This demonstrates that the limit 5i(^) — > | is quite continuous when viewed at finite times with well-defined 
physical initial conditions, although the late time limit is discontinuous. 
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FIG. 4. From top to bottom the curves are the trace, pressure, and energy density for m = and f = 0. The state is a fourth 
order adiabatic state at t = 1 with n\ — — 2 and nk — for all k > 2. The trace, pressure, and energy density approach 
their AF fixed point values at late times. 
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In the case v = | there are two distinct behaviors depending on whether m and £ are separately vanishing or not. 
In the massless, minimally coupled case we proved in Section IV that the Allen-Folacci value is a fixed point at late 
times. In Figure 4 we show the approach of the energy density, pressure, and trace to their Allen-Folacci values for 
the massless minimally coupled field. The field is in an "n-particle" state with n\ = ?i2 = 2 and — for all k > 2. 
In Figure 5 we show the behavior of the energy density and trace for the case m — 0.3 and £ = —0.0075 when the 
field is in a fourth order adiabatic vacuum state. 
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FIG. 5. From top to bottom the curves are the trace, pressure, and energy density for m = 0.3 and £ = —0.0075, = and 
the fourth order adiabatic state at the initial time t = 1. The curves behave linearly in comoving time t at late times, with 



their slopes given by Eqs. (4.29a) and (4.29t) 
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FIG. 6. From top to bottom the curves are the trace, pressure, and energy density for m = 0.1 and £ = —0.3, nt — and 
the fourth order adiabatic state at the initial time t = 1. 



When > | the analysis at the end of Section IV shows that the leading order state dependent terms in the 
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energy-momentum tensor will generally grow exponentially with time. This is illustrated in Fig. 6. The exponential 
growth that occurs here is similar to the well known one for the classical scalar field when m and £ have values such 
that v > | |J8|,|9) . Such fields are tachyonic and presumbably of little physical interest unless interactions are added 
to stabilize them. 



VI. INFRARED SCALING AND THE GENERALIZED ANOMALY 

We have found that for all 5i(z/) < § the renormalized expectation value of (T a b) approaches the de Sitter invariant 
Bunch-Davies value for any physically admissable initial state, whereas it approaches the de Sitter invariant Allen- 
Folacci value for any physically admissable initial state in the massless, minimally coupled case. Since all the initial 
state dependence vanishes asymptotically, these state independent de Sitter invariant fixed point values for {T a b)ren 
must be purely geometrical in origin. Indeed, both the Bunch-Davies point-splitting calculation and the Hadamard 
calculation of Allen-Folacci rely only on the properties of the two-point function of the scalar field G(x, x') for x' — > x. 
Hence, the BD and AF asymptotic values of (T a b)ren are certainly "pseudo-local" in the terminology of Ref. jfO|, i.e. 
they are expressible in terms of purely local functions of the RW scale factor a(rj) and its derivatives. 

If we specialize now to zero mass, m = 0, then on simple dimensional grounds the asymptotic (T a b) r en can be 
expressed purely in terms of local conserved tensors of fourth adiabatic order. Although we have used adiabatic 
subtraction methods to renormalize (T a b) it is known that the value of (T a b) ren so obtained is equal to that in a 
fully covariant procedure such as dimensional regularization or covariant point-splitting p7| , p8| . In a fully covariant 
procedure, which yields a local conserved tensor of fourth adiabatic order, only the Riemann tensor together with its 
covariant derivatives and contractions can appear. Hence (T a b) r en for m = must be expressible entirely in terms of 
such local geometrical tensors. 

In four dimensions the only such local tensors are linear combinations of ^'H a b, ^ H a b, and ^H a b, where 



(i) 



2g ab UR - 2V a V b i? + 2RR ab - ^g ab R 2 , (6.1) 



vanishes in de Sitter spacetime and 



2 1 1 

^'H a b — R c a R C b — -^RRab — ^RcdR cd gab + -^R 2 g a b ■ (6.2) 

In RW spacetimes, which are all conformally flat, the tensor ^H a b is proportional to V->H a b and hence vanishes as 
well. Therefore the only non- trivial fourth order conserved geometrical tensor in de Sitter spacetime is ( 3 >H a b and we 
conclude that the fixed point BD and AF values found in our previous analysis are proportional to 

R 2 „ f R x 2 

V " ; Hab — 

Furthermore, since 



W Hab = ^g ab = 3(^) g ab . (6.3) 



{3) H ab g ab = -R ab R ab + ^ = \ (RabcdR abcd - m ab R ab + R 2 ) = j, (6.4) 

in RW spacetimes (where the Weyl tensor vanishes), the coefficient of H a b is proportional to the coefficient of the 
Gauss-Bonnet integrand in the trace of (T a b). Such a term in the trace is known to correspond to a non-local but 
nevertheless fully covariant action and this action is precisely the same as that generated by the trace anomaly of free 
conformal fields |po| . Since we have obtained fixed point results for the asymptotic values of {T a b) ren for massless 
fields in de Sitter space which are purely geometrical and proportional to ^H a b, even for non-conformal massless 
fields, we can give a definite meaning to the value of the proportionality coefficient and the non-local anomaly-like 
term in the effective action even when £ ^ i. 

Let us define the generalized anomaly coefficient by fixing the normalization 

Urn (T ab ) ren = -^- 2 ^H ab = (£\ g ab . (6.5) 



21 



With this normalization we find from the asymptotic value of {T a b) ren for a scalar field in de Sitter space that 
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(6.6) 



The value of Q 2 for a conformally invariant field (m = 0, £ = i) is ^gg anc ^ corresponds to the pure trace anomaly 
coefficient. The first member of (3.6) provides the generalization of this coefficient away from the conformal case. The 
discontinuous behavior at £ = has been discussed in Ref. . As we have seen it arises from the singular behavior 
of the spatially constant zero mode of the massless, minimally coupled field, which is non-oscillatory and hence cannot 
be quantized as a Fock mode in the same fashion as the oscillatory modes. We discussed this discontinuity in detail 
in Section IV. 

The connection of the tensor ^H a b with the trace anomaly may be seen from the general form of the effective 



action for the anomaly in a conformally flat space with metric 20 



9ab 



e 2a g a b , 



(6.7) 



namely 



1 (r* 



crA i a+- I G- - .09.09 R ] a 



(6.8) 



\RU 2 + \{V a R)V a is the unique fourth order differential o pera tor acting on scalars 



8) can be obtained by 



where A 4 = d + 2i? Qb V a V 6 

which is conformally covariant. A fully covariant but non-local form of the effective action 

solving \J—g (G — §LJ-R) = y/— 5 (G — §D-R) +4-^/— gA.4<r for a. In that non-local form all reference to the separation 
of the metric into background and conformal factor as in ( |6 . 7| ) disappears. 



The energy-momentum tensor following from the variation of the local form of the effective action (6 
to the background metric g a b, is given by Eq. (2.9) of Ref. 
considerably on the Einstein space R x S 3 , where G 



with respect 



] . The form of this energy- momentum tensor simplifies 
09.09 R = 0. If we set a = log a(rf) then this is equivalent to 



evaluating the components of this energy-momentum tensor in a general RW spacetime with closed spatial sections 
(k = +1). Using the expressions for W H ao and ^H a t, in a general RW space in terms of a(rf) and its derivatives, we 
quickly find that 
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(6.9) 



Hence the tensor ^H a j, which is called "accidentally conserved" in Ref. [Q 
a non-local covariant effective a ctio n related to the trace anomaly, which has the local form (6 
conformally decomp osed as in (3.7 
The last term in (| 



is associated in fact with the existence of 
when the metric is 



■9D would not have been present had we varied the fully covariant but non-local form of the 
anomal ous effective action, in which all reference to the background g a b drops out. Equivalently, it is just canceled if we 
add to (6.9) the Casimir energy on R x S 3 , which is determined by the anomaly by a further conformal transformation 
from flat space jOl. In either case, we should drop this last term in (3.9) which depends on the arbitrary background 
g a b- Evaluating the first two terms on de Sitter space we find that it is exactly the value of the asymptotic form of the 
renormalized energy- mo ment um tensor (T a b) in de Sitter spa ce (| 6.5[ ), found previously, with the value of Q 2 for the 



massless scalar given by (|6.6|). Thus, the effective action (3.8), associated with the conformal trace anomaly, appears 



in the effective action for a massless scalar field, even for non-conformal couplings, £ ^ 

Since it is associated with the anomaly, the physical significance of S g in the quantum effective action for the 
scalar field is that it determines the scaling behavior of the field theory under global Weyl transformations of the 
background space. Using the fact that the Euclidean effective action is given by I e g = —S c g and that the Euclidean 
continuation of de Sitter space is S 4 with Euler number x = 2, we can vary I c g with respect to a constant a = ao 
and obtain 



dies _ digg _ Q 
dao da 



£-2 [ d 4 x^G = Q 2 X = 2Q 2 , 

o27T z J S i 



(6.10) 



dcrn 



dl p gee is a global rescaling of the S 4 . As a check, this relation can be verified explicitly by the C function 



evaluation of the Euclidean effective action, 



1 



/ eff = - Tr log I 



.09.09 



2ds 



(6.11) 
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where the generalized zeta function for the Euclidean continuation of the wave operator appearing in the Tr log is 



COO = V d n \~ s , (6.12) 



ri=0 



in terms of its eigenvalues X n with degeneracy d n on S 4 . This sum is convergent for -Ji(s) > 2 and defines a meromorphic 
function of s which is analytic near s = 0, where its derivative is required. Introducing a mass scale \x to keep £(s) 
dimensionless for all s leads in a standard calculation to 11421 



1^ 

2ds 



= -C(0)log(/za)+Ii(*/), (6.13) 



where /i(^) is a certain finite function of v, which from the definition of v (3^) becomes independent of a when 
to = 0, and the value £(0) is given by |^] 

-«w-5(-*' + i« , + ffii)-B-5«- 1 »'- ( " 4 » 



when 77i = 0. By making use of Q, ( |6~Tl| ), and ( |6~T3| ), we find that 



a^-J eff (m = 0) = -C(0) = 2Q 2 , (6.15) 
da 

is the behavior of the ef fectiv e action for a massless scalar field with £ > under global Weyl rescaling of the metric, 
exactly as predicted by( |6.10 ) and the previous discussion based on the anomalous action (6.8). 



When to = and £ — > 0, the integral representation of the function 1\ ( v) de velops a logarithmic singularity, which 
can be traced to the vanishing of the n = eigenvalue in the expression |Tl|) for £(s). In this case the n = mode 
must be excluded from the ultraviolet regulated sum over modes, which has the effect of adding one unit to £(0) in 
the infrared scaling behavior of the effective action p3| ], and accounts for the addition of — | in Q 2 in the minimally 
coupled case. Hence the discontinuous behavior of Q 2 found in ([Tf]) by our analysis of the asymptotic attractor 
behavior of the energy-momentum tensor in de Sitter space is precisely the same as that occuring in the effective 
action under global Weyl rescalings when to = and £ — > 0. 

Note also that the global Weyl variation is given in terms of the trace of the energy-momentum tensor by 

a—I eS (m = 0) = -[ d 4 x Jg (T) ren = -^-a 4 (T) ren . (6.16) 
da J S 4 S 

Since the £ function method is fully covariant, the renormalized trace (T) ren computed in this way cannot contain 
non-covariant contributions, and must be expressible entirely in terms of local curvature invariants. 

Finally we may consider relaxing the condition to = 0. If m 2 > 0, the asymptotic form of the ener gy- momentum 
tensor for an arbitrary initial state is given by the BD value. However, if we expand the BD result (|3.6|) in powers 
of R/m 2 , we find that it contains no adiabatic order four R 2 terms, beginning instead with R 3 /m 2 . Mathematically 
this is because all terms up to fourth adiabatic order have been removed by the ultraviolet regulating procedure of 
point splitting or adiabatic subtraction. Hence the coefficient of ^H a b at asymptotically late times in an arbitrary 
physical state is given by 

Q 2 = for to 2 > , (6.17) 
and no anomalous S c g term appears in the quantum effective action for a massive field. This is consistent with the 



interpretation of the anomalous term (6^) in the effective action of the scalar field as an infrared effect, since the 
fluctuations of a massive field decouple at large distances or late times, and should induce only strictly irrelevant 
operators in the effective action in the far infrared, which are suppressed by positive powers of R/m 2 . 

Only in the strictly conformal case, m = and £ = g is the infrared effective action equal to that obtained by 
ultraviolet methods, such as the ci2 coefficient in the Schwinger-DeWitt proper time expansion. However, our analysis 
of the fixed point behavior of (T a b) rerl in de Sitter space shows that its coefficient is connected with the global or 
extreme infrared scaling of the effective action, and this asymptotic behavior does not depend on the field being 
conformally invariant. The asymptotic attractor behavior of the energy- momentum tensor in de Sitter space defines 
an infrared scaling coefficient that reduces to the trace anomaly coefficient in the conformal case, but is a much more 
general concept than the trace anomaly coefficient, since it is well-defined for all massless fields, conformal or not. It 
is well-defined even for massive fields, although as ( |6. 17 ) shows, it vanishes in this case. 
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VII. DISCUSSION AND CONCLUSIONS 



In this paper we have considered a quantum scalar field in a fixed de Sitter background. We have studied the 
late time behavior of the renormalized energy-momentum tensor and have found two important cases in which, for 
arbitrary physically admissable states, the energy-momentum tensor approaches a particular value at late times. The 
values approached are those of the energy-momentum tensor in the Bunch-Davies and Allen-Folacci states. Thus, in 
this sense, these special quantum states behave as fixed point attractors. 

In the case < § we have shown that for all fourth order adiabatic states that are infrared finite the energy- 

momentum tensor approaches the BD value at late times. The longest time scale for the state dependent terms to 
redshift away is t = (3 — 2v)~ 1 . This has been numerically verified for various values of m and £ when the fields are in 
a fourth order adiabatic state. For the case in which < | — v «1 and m and £ are small, we numerically observe a 
more complicated behavior. The energy-momentum tensor quickly approaches the AF value and grows then linearly 
with comoving time. Our analytic proof implies that it must then slowly decay to the BD value. It is worth noting 
that for Si(iv) > the redshift of the state dependent terms in the quantum expectation value of T a b is slower than 
one might guess from the redshifting of classical matter or radiation, i.e. a~ 3 or a~ 4 , respectively. 

For the case v = | we have to distinguish two different possibilities. If the field is massless and minimally coupled 
then we have proven that, for all fourth order adiabatic states that are infrared finite, the energy-momentum tensor 
approaches the AF value at late times. This is true for both vacuum and initially populated states. For any other 
values of m and £ the renormalized energy-momentum tensor grows linearly with comoving time, indicating that 
back-reaction effects need to be taken into account. The sign of the linear growth depends on the sign of £. 

For the tachyonic cases v > | there is no attractor state. Instead, for most values of m 2 and £ the renormalized 
energy-momentum tensor grows like a 2iy_3 at late times for all physically admissable states. Thus back-reaction effects 
again need to be taken into account. 

For the cases in which either the Bunch-Davies or Allen-Folacci states serve as attractors in the above sense and the 
mass, m of the field is zero, we have shown how these results are connected to the appearance of a certain non-local 
term in the quantum effective action for the scalar field. This term gives rise to the local geometrical tensor ^H a b 
in the asymptotic form of (T a b) at late times, and also determines the global scaling behavior of the effective action 
for massless fields. Determining this scaling behavior and relating it to the asymptotic (T a b) in de Sitter space has 
allowed us to generalize the notion of the trace anomaly to massless, non-conformally coupled scalar fields, in the 
sense that the coefficient of this non-local term in the effective action is well defined even for £ ^ | . 

The interplay between the UV and IR properties of the state and the mode sums contributing to the energy- 
momentum tensor is a theme running through all of these results. We had to be careful to remove the UV divergences 
from the unrenormalized (T a b) in order to analyze the late time limit. Since we have found state independent de 
Sitter invariant results in both the v < | and massless, minimally coupled cases, and since all the state dependence 
resides in the finite k modes, their form at high k being restricted by the requirement of matching the adiabatic order 
four vacuum, it is clear that state independent results for (T a b) r en are possible only because the contribution to the 
BD or AF expectation value comes from arbitrarily large k at very late times. In fact, inspection of the renormalized 
expectation value of (T a b) expressed as a mode sum, after the fourth order adiabatic subtraction has been made, 
shows that the finite contribution comes from k ~ a, i.e. when the physical wavelength of the mode is of order of the 
de Sitter horizon. At very late times, this corresponds to arbitrarily large values of the coordinate wave number k. 

The finite difference between the BD value and the AF value in the m = £ = case comes entirely from the k = 1 
mode in closed spatial sections, which is a purely IR effect. This leads to a finite discontinuity in the infrared scaling 
properties of a massless field since the value of the energy-momentum tensor at m = £ = is different from its value 
in the limit m = and £ — > 0. The appearance of the ( 3 'H a h tensor and the corresponding non-local action is quite 
generic for massless fields; its coefficient Q 2 vanishes only if the mass is non-zero. Hence we should expect that, 
although they are certainly not conformal, gravitons will also contribute to this same infrared effective action with a 
finite value of Q 2 , which can be determined in the same way by a background de Sitter calculation of their quantum 
{T a b) at late times. We plan to present the results of this calculation in a future publication. 

It is also interesting to note that the coefficient of the generalized trace anomaly, Q 2 , is not generically positive, in 
contrast to all the previously known examples of massless conformal fields |E3 . It appears that the reason for this 
is that a positive Q 2 comes from the ultraviolet behavior of (T a b) for conformal fields, while the infrared behavior of 
non-conformal fields can contribute a negative value. For any Q 2 ^ the new term in the effective action leads to 
dynamics for the RW scale factor which is quite different from the Einstein theory, and remains to be investigated in 
a full dynamical back-reaction calculation. 
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APPENDIX A: THE CASES OF IMAGINARY AND INTEGER VALUES OF v 



In this appendix we show that the proof in Section III works for imaginary values of v as well as for v = 0, 1 . For 
i magin ary values of v it is useful to write v = ij with 7 a positive real number. Then the formulas in Eq. ( 3.23a ) and 
( |3.23b| ) are the same as before. However the values of the Ai are different. They are given by 



At = 

A 2 
A 3 



2k 



coth(77r)csch(77r)((l + 2n/c)|c2| 2 + n k ) - -coth(77r)csch(77r)(l + 2n k )(cic 2 + c*c 2 ) 



--csch(77r)(l + 2n k )(cic* 2 - c*c 2 ) 



~k 6 



2k 



(1 + csch 2 ( 7 7r)((l + 2n k )\c 2 \ 2 + n k ) - -csch 2 ( 7 7r)(l + 2n k )( Cl c* 2 + c*c 2 ) 



coth(77r)csch(77r)((l + 2n k )\c2\ 2 + n k ) - -coth(77r)csch(77r)(l + 2n k )(cic 2 + 0*02) 



+ -csch(77r)(l + 2n k ){cic* 2 - c\c 2 ) 



(Al) 



The last terms in the expressions for A\ and A3 are imaginary. Substitution into Eqs. (3.23a) and (3.23b) shows that 
the contributions of these terms to the energy-momentum tensor cancel. 

Noting that v is imaginary, one sees that the argument for the vanishing of the first integral in Eq. ( 3.27 ) is 
unchanged. Since the terms in the integ rand are b eing bounded for the second integral one must take the real part 
of j3i and substitute that for (3i in Eqs. ( 3.28| ) and (3.31). After making this substitution it is still the case that each 
term in ( 3.3 ip vanishes in the limit rj — > 0~. The argument for the vanishing of the third integral in Eq. fl3.27 ) is 
unchanged when v is imaginary. Thus, the proof is valid for imaginary values of v. 

For v = n = 0, 1 we use the identity 

H£\z) = J n (z) + iN n (z) , 

Hi 2 \z) = J n {z) - iN n (z) , (A2) 
and then use the well known power series solutions for Nq and N± to write 



N n {z) = -J n {z) log(z) + P n {z) , 

7T 



where P n (z) is a series of the form 



Pn{z) = 

Because of the log(z) term it is useful to write 

dN n (z) _ 2 dJ n (z) 



2j-n 



dz 



with Q n {z) a series of the form 



7r dz 



.(*) = !/ 



log(z) + Q n (z) , 



y 2j-n-l 



(A3) 



(A4) 



(A5) 



(A6) 



We can conclude that Eqs. (3.23a) and (3.23b) remain the same but with different expressions for A4, Si, and fy. The 
new expressions for the Ai are given by 



A, =~ 



7T 

2k: 



l + ^(log Z y ((l+2n fc )|car+nfc) + - 1- ^(logz)^ (1 + 2n k )( Cl c* 2 + c\c 2 ) 



2i 

— log z (1 + 2n k ){c x c* 2 - c*c 2 ) 



A 2 
A 3 



IT 

k 

7T 

2k 



log2;((l + 2n fc )|c 2 | 2 +n fe ) 



1 



1 i 

- logz (1 + 2n k ){dc* 2 + c*c 2 ) + -(1 + 2n k ){ac 2 - c\c 2 ) 

7T A 



((1 + 2n k )\c 2 \ 2 + n k ) - -(1 + 2n k )( Cl c* 2 + c*c 2 ) 



(A7) 



2G 



Note that the Ai are now functions of both k and log(z). The new expressions for the fti and Si are given in Table 2. 



i 


ft 


Si 


1 


5 + 2n 


z 5 J 2 (z) 


2 


5 


z 5 J n (z)P n {z) 


3 


5 - 2n 




4 


3 + 2n 


zMZ(z) 


5 


3 


z 3 J n (z)P n {z) 


6 


3 - 2n 


z*PZ(z) 


7 


3 + 2n 




8 


3 


z 4 (£j n (z)) P n (z) + J n (z)Q n (z)) 


9 


3 - 2?i 


2z i P n (z)Q n (z) 


10 


3 + 2n 




11 


3 


z 5 (£j n (z))Q n (z) 


12 


3 - 2?i 





Table 2 



Substitution these expressions into Eqs. ( 3.23a ) and ( 3.23b ) shows that for the first integral in ( 3.27 ) the expressions 
are of the same form, except that some term s have factors of log(z). However these do not prevent the terms in the 
first integral on the right hand side of ( |3.27 ) from vanishing asymptotically. For the second integral there are still 
terms of the form given in Eq. fl3~2§| ). However, some of them also have factors of log(z). Inserting factors of log(z) 
into Eq. (3.31) and computing the integrals, one finds that the terms all vanish in the limi t t] — > CP. The factors of 
log(z) also do not affect the asymptotic vanishing of the third integral on the right in Eq. ( 3.2 7| ) . 
Therefore, in all cases where Si(z/) < | the quantity (T a b)sD vanishes in the limit r; — ► 0~. 



APPENDIX B: THE HARMONIC OSCILLATOR WITH u -> 



In this Appendix we give a detailed discussion of the harmonic oscillator with vanishing frequency, pointing out the 

3 
2 



analogy with the v — > | limit in de Sitter space. Consider a simple harmonic oscillator with Hamiltonian 



H = \tf + i^ 2 . (Bl) 
The single degree of freedom <f> can be quantized by introducing the operator representation 

<t>(t) = (ct>(t))+am + aH*(t), (B2) 
where and a are creation and destruction operators obeying 

[0,0+1 = 1. (B3) 

The cano nical commutation relations [</>, cf>] = i are satisfied provided the mode function i\) obeys the Wronskian 
condition (2.S). The equation of motion for the Hciscnberg operator <f>(t) implies that the mode functions also satisfy 
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4> + u; 2 ip = 0. (B4) 

In order to satisfy the Wronskian condition we choose the fundamental normalized complex solution to this equation 
to be 

f(t) = -±= e-^ . (B5) 



The general solution satisfying both the Wronskian condition and (B4) can be written as the linear combination 

= <*/(*) +/3/*(t), (B6) 

where 

H 2 -|/?| 2 = l. (B7) 

Thus, up to an irrelevant overall phase, the Bogoliubov coefficients can be parameterized in terms of two real param- 
eters in the form 

a = cosh 9 , 

/3 = sinh0 e~ lS . (B8) 

We could keep the c-number expectation value (4>{t)) non-zero in general, which would lead to the general Gaussian 
state with displaced origin that we considered in some earlier papers Since we do not require it in our free field 
de Sitter calculation, we specialize to the case where (4>(t)) = (4>(t)) = in the following. With this restriction the 
field operator can be written simply as 

4>(t) = ai>(t) + aty*(t). (B9) 
We make use of the freedom in the a and (3 coefficients to require 

(aa) = (aV) = , 

(aa f ) = (a f a) + 1 = , (BIO) 

with no loss of generality. The positive real parameter a > 1 has the interpretation a = 1 + 2n, with n the average 
number of "particles" in the a basis. Thus we need the three parameters 9, 5, and a to specify the general semi-classical 
(coherent) state of the field. 

Using the definitions above it now follows that the average energy in this state is 

E = <ff) = |a(l + 2|/3| 2 ). (Bll) 

Note that we have not insisted that the state be an eigenstate of the number operator a^a so n is only the average 
number of particles in the state, which can take on any non-negative value and need not be an integer. If n > and 
the state is a Gaussian, consistent with all connected correlations vanishing except the two-point function, then it 
turns out that the state is necessarily a mixed state || . 

We are interested in the nature of the "vacuum" state in the limit lu —> 0, i.e. when our single degree of freedom 
becomes that of a free particle. In this limit we can no lon ger r etain the complex oscillating mode function basis since 
their normalization diverges due to the (2oj)~2 factor in (|B5|). However, the mode equation if) = clearly possesses 
the regular real solutions u = 1 and v = t, which we can obtain from the complex oscillatory solutions by taking 
appropriate linear combinations in the limit of vanishing frequency, i.e. 

^ {Vf(/ + /•)} = !, 
v(t) = i hm(-L(/ -/*))=t. (B12) 



o I V2cj 



These definitions are analogous to Eq. (4.7) of the text, with u> replacing 3 — 2^ 



The general linear combination of modes can then be rewritten in the form 
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'2uj 



i—(a - 0)v = Bu + Av -> At + B , 



(B13) 



where the quantities 



B = lim 



1 



are analogous to those defined by (4.8) of the text. They also satisfy 

A*B-B*A = i(\a\ 2 -\0\ 2 )=i. 
We can also define the time-independent Hermitian operators 



(B14) 



(B15) 



1 



[(a +0) a+(a + (3)* a f ] -> Ba + B*a) 



Q 



P = -i x j- [(a -p)a-(a- 0)* a 1 ] -> Aa + A*a f 



analogous to ( 4.11 ) in the de Sitter case. They obey the canonical commutation relations 

[Q,P]=i, 

so that in the limit lo — > 



and in the same limit 



E 



<f>{t) -» u(t)Q + v(t)P = Q + tP, 



; <P 2 ) = 5lim> \a-(3\ 2 } = 



(B16) 



(B17) 



(B18) 



(B19) 



Since a is just an overall factor and we are interested in pure vacuum-like states, we can set a = 1. 

We see that depending on exactly how we take the limit lo — > 0, we can get many different values for the energy. 
If we tried to keep the pure positive frequency solution tp = f, i.e. a = 1 and (3 — 0, then the energy would vanish 
in the limit lo — > 0. However, by computing the correlation function, (<f)(t)<j){t!)), we would find that the correlator 
diverges in this limit. In our simple example it is clear why. The ordinary ground state vacuum wave function of the 
simple harmonic oscillator is a Gaussian with a width proportional to Hence (Q 2 ) ~ lu^ 1 in this state, which 

becomes non-normalizable in the limit, unless we were to put the particle in a box of finite volume. But this option 
is unavailable to us if the "particle" is a mode of a field and the range of the field variable <f> is (— oo,oo). Thus, we 
must reject the positive frequency harmonic oscillator vacuum in the limit lo — > since it is a non-normalizable state. 

In order for the state to remain normalizable and the mode func tions bounded, the quantities A and B must 
remain finite in the limit of vanishing lo. Inspection of the definitions (B14) shows that this requires that the original 
Bogoliubov coefficients \a\ and \/3\ go to infinity. In fact, using (BS) in the definitions of A and B a simple calculation 
shows that the Bogoliubov parameters and S must behave like 



-iiog(f)+iog|4 

7T + LOT , 



(B20) 



as lo — > 0, so that 



A = -i\A\, 
B 



and 



2jA\ + 2-^ T > 



(B2I) 



remain finite in that limit. Thus, the two complex numbers A and B obeying (B15) and the state they characterize 
can be specified by the two real numbers \A\ > and T, up to an unobservable overall phase. Comparison of these 
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expressions with the corresponding equations (4.7) and (4.13) in the paper by Allen and Folacci |l7| shows that our 
A here is equivalent to Allen and Folacci's 2A. This is because of a factor of 2 difference in the normalization of their 
time dependent mode v(t). 

Since 9 — ► oo, in a sense the normalizable "vacuum" state with finite (P 2 ) = \A\ 2 and finite (Q 2 ) = \B\ 2 is infinitely 
far from the usual ground state vacuum of the harmonic oscillator, and is not an eigenstate of the free particle 
Hamiltonian at u> = 0. Instead, it is just a Gaussian state centered at (Q) = and (P) = with normalized wave 
function || 

1 \ A \ ..2 



* (4) = (i^ ex H-^' J J' < B22 » 

in the position representation where Q^f — and = — i^5 f . This state is time reversal invariant if and only 
if T = 0, in which case B — becomes real. On the other hand, an eigenstate of the free Hamiltonian P 2 jl is 
necessarily non-normalizable since diagonalizing P requires infinite spread in Q. These remarks carry over equally 
well to the k — 1 mode in the de Sitter case. 

There is no reason to diagonalizc H in the field theory case as Allen and Folacci do in their Eqs. (4.26) to 
(4.28) |l7j ], since the energy-momentum tensor is time dependent in general, and back-reaction has been neglected 
in our simple calculations. Diagonalizing just the matter Hamiltonian without taking into account the metric is no 
better than a leading order mean field approximation to the much more difficult full problem with metric fluctuations 
taken into account. Given that what we are doing is only a semi-classical approximation to this full problem in any 
case, what makes the most sense from our perspective is to take well-defined normalizable initial adiabatic states for 



the scalar field and allow them to evolve as they will. Our initial state trial wave functional (B22) then remains a 
bounded Gaussian functional, consistent with semi-classical mean field methods (although it may well spread over 
time depending on the dynamics). 

As should be clear, the non-normalizable positive frequency vacuum with a = 1 and (3 = is analogous to 
the Bunch-Davies vacuum, while the normalizable state related to it by an infinite Bogoliubov transformation (but 
only a finite shift in zero point energy) is analogous to the Allen- Folacci states parameterized by A and B or \A\ 
and T. However, wher eas the energy of the normalizable state is finite and arbitrary, depending on the arbitrary 



\A\ 2 coefficient in (B19), and the non-normalizable positive frequency state has zero energy in this simple harmonic 



oscillator example, the converse is true at late times in de Sitter space. This is due to the different kin emat ics and 



form of the energy-momentum tensor of the scalar field in the de Sitter case. Comparison of the energy ( B19 ) of the 
harmonic oscillator example and the corresponding k = 1 late time energy density in the de Sitter case with m = 
and £ small, viz. 

shows that the difference is the redshift factor of a~ 2 multiplying the infrared finite {P 2 } in the de Sitter case. Thus, 
whereas this can take on the finite state dependent value \A\ 2 in the harmonic oscillator example, its asymptotic late 
time value vanishes in de Sitter space when multiplied by a~ 2 , and becomes state independent in any infrared finite 
state. Further, in the simple harmonic oscillator the usual ground state becomes an eigenstate of P and has zero 
energy as u> — > 0, since lu 2 ((Q + vP) 2 ) ~ u — > drops out of the energy in this limit, whereas in the de Sitter case, 
since £ goes to zero only linearly with 3 — 2v, the Bunch-Davies expectation value receives a finite contribution from 



the k = 1 mode from the infinite spread in £(Q ) which has no a redshift factor in (B23) 
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